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Abstract 



We compute the pion inclusive momentum distribution in heavy-ion col- 
lisions at AGS energies, assuming thermal equilibrium and accounting for 
density and expansion effects at the time of decoupling. We compare to data 
on mid rapidity charged pions produced in central Au + Au collisions and 
find a very good agreement. The shape of the distribution at low m; - m is 
explained in part as an effect arising from the high mean pion density achieved 
in these reactions. The difference between the positive and negative pion dis- 
tributions in the same region is attributed in part to the different average 
yields of each kind of charged pions. 

PACS number: 25.75.-q 
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A great deal of experimental effort has been devoted in recent years to the production of 
a highly compressed state of matter in high-energy collisions between heavy ions. It has been 
speculated that at sufficiently high baryonic densities or temperatures, a phase transition 
from hadronic matter to a quark-gluon plasma can occur. The study of particle spectra 
should provide useful information about the dynamics and evolution of the kind of matter 
formed in these reactions. 

Some of the most copiously produced particles in central high-energy heavy ion colli- 
sions are pions. An understanding of pion production during these collisions has long been 
searched, particularly in view of one of the most remarkable properties exhibited by their 
spectra commonly referred to as an enhancement of either the low or the high transverse 
momentum regions in the inclusive single pion distribution, as compared to p-p collisions. 
This property is concomitant with the difficulty to describe the invariant transverse mass 
distribution with a single exponential function ||. 

Several possibilities have been put forward to explain the peculiar shape of the distribu- 
tions, amongst them, the different contributions to the pion yield coming from the decay of 
A resonances produced during different stages of the collision [[| and the superposition of 
primary pions and pions coming from resonance decay, mainly A's [|J. The importance of 
transverse flow in describing the spectra has also been stressed j5|. The proper treatment 
of Coulomb final state interactions has also been pointed out. More recently, Wong and 
Mostafa noticed that when particles feel the effects of a boundary during the evolution 
from the first stages of the collision to their final free streaming, their momentum distribu- 
tion is affected due to the discretization of energy levels and the correspondingly different 
density of states introduced by the finite size of the system just before freeze out. This 
idea is based on the concept of a pion liquid, first discussed by Shuryak |7| and was further 
developed in Ref. || with the introduction of a finite chemical potential associated to the 
mean pion multiplicity per event in central collisions. 

Calculations based on the effects of a boundary resort to the assumption of thermal 
equilibration. They are successful in reproducing the concave shape of the distribution at 
high transverse momentum but fail to describe its overall fall off. This failure could have 
been anticipated since, as pointed out in Ref. 0, thermodynamics is tantamount of hydro- 
dynamics and a simple thermal spectrum should be corrected by the Doppler shift resulting 
from collective expansion. A challenge remained as to how to incorporate hydrodynamical 
effects into a description of the transverse momentum distribution based on discrete energy 
states. In this work we meet this challenge and lay down the formalism required to account 
simultaneously for expansion and boundary effects in a phenomenological calculation of pion 
spectra in relativistic heavy-ion collisions. We apply this formalism to compute the invariant 
transverse momentum distribution for mid rapidity pions at AGS energies. We find a very 
good agreement with data on central Au+Au reactions at 11.6 A GeV/c. Further distortions 
of the spectra arising from Coulomb effects will be discussed in an upcoming work. 

In order to collect the necessary ingredients for the calculation, recall that if the average 
pion separation d ever becomes smaller than the average range of the pion strong interaction 
d s (~ 1.4 fm) during the evolution of the collision, the pion dispersion curve can be modified 
and the collective properties of the pion system resemble those of a liquid rather than those 
of a gas |7j. An important consequence is the appearance of a surface tension that acts 
as a reflecting boundary. The E-802/866 collaboration has reported that a baryon density 
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of about eight times normal nuclear density is achieved in central Au+Au collisions at 
11.6 A GeV/c. A large fraction of this density is due to pions. From Fig. 3 of Ref. f2|, one 
can read that the total number of charged pions one unit around central rapidity in this kind 
of reactions is about 200. Under the assumption that the number of neutral pions in the 
same rapidity interval is half the total number of charged ones, the total pion yield in central 
collisions at mid rapidity is dN^/dy ~ 300. Since the average pion separation is inversely 
proportional to one third of the pion density n n = (l/Ato)dN n /dy, where A(~ 64 fm 2 ) is the 
transverse area of the reaction and to(~ 1 f m ) is a typical formation time fTO"! , d ~ 0.6 fm 
< d s and thus the condition to regard the pion system as a liquid is met. 

Pions that move towards the boundary of the system feel the attractive potential behind 
them and are reflected back. The reflection details depend on the wave length of the given 
particle but the property introduced by the reflecting surface is that it allows very little 
wave function leakage and to a good approximation, the pion wave functions vanish outside 
the boundary. When, as a consequence of the expansion of the initially compressed and hot 
system, the pion average separation becomes larger than the range of strong interactions, 
the system becomes a free gas but the transition between the liquid and the gas phases is 
very rapid and the pion momentum distribution should be determined by the distribution 
just before freeze out. 

Thus the system of pions can be considered as confined and their wave functions as 
satisfying a given condition at the boundary just before freeze out. In this case, the energy 
states form a discrete set. The shape of the volume within the confining boundary deserves 
some attention. For reactions with a large degree of transparency, Bjorken like geometry, 
with a predominantly longitudinal elongation, seems better suited. However, for AGS ener- 
gies, a significant amount of stopping has been reported by the E-802/866 collaboration in 
central Au+Au reactions 0. In this case, a more symmetric geometry between transverse 
and longitudinal directions seems appropriate. We thus consider a scenario in which the 
system of pions of a given species is in thermal equilibrium and is confined within a sphere 
of radius R (fireball) as viewed from the center of mass of the colliding nuclei at the time 



of decoupling. As discussed in Ref. [11], this time needs not be the same over the entire 



reaction volume. Nevertheless, in the spirit of the fireball model of Ref. |LIJ, we consider 
that decoupling takes place over a constant time surface in space-time. This assumption 
should be essentially correct if the freeze out interval is short compared to the system's life 
time. 

In the absence of expansion, the solution has been found in Ref. ||. This involves solving 
the Klein-Gordon equation to find the stationary wave functions satisfying 



d 2 



V(r,«) = 0, (1) 



subject to the condition 



1>(\r\=R,t) = 0, (2) 



and finite at the origin. In order to incorporate the effects of a hydrodynamical flow, we 
observe that the presence of an ordered motion, represented by a four-velocity field ii M = 
7(r)(l, v(r)), amounts to a redistribution of momentum in each of the fluid cells, as viewed 
from a given reference frame (the center of mass in our case). The tendency of matter to 
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occupy a larger volume is compensated by the distribution of momenta in each cell becoming 
narrower 



|12| . The distribution in the cell becomes also centered around the momentum 
associated with the velocity of the fluid element. Consequently, the thermal spectrum in 
each cell should be described on top of this collective flow, that is, referred from the collective 
fluid's element momentum. To describe this behavior of the expanding, bound system of 
pions, we make the substitution of the momentum operator p M by — mu^, where m is 
the pion mass. The term mu^ represents the collective momentum of the given pion fluid 
element. The corresponding equation becomes 



d , „\2 



i— 7717(7")) + ( — z'V — m/y(r)v(r)) + m 2 



1>{r,t) = (3) 



and we look for stationary solutions subject to the same condition in Eq. (fj) and also finite 
at the origin. We consider a parametrization of the three velocity vector v(r) that scales 
with the distance from the center of the fireball. 

v(r)=/^f, (4) 

We identify this velocity with the transverse flow velocity and ignore any asymmetry between 
transverse and longitudinal expansion. < 3 < 1, represents the surface fireball velocity. 
The corresponding explicit expression for 7(7") is 

7(r)= , 1 (5) 
" -3 2 r 2 /R 2 

Eq. (0) with the gamma factor given by Eq. (|5J) can only be solved numerically. In order 
to provide an analytical solution, we approximate the function 7 by the first terms of its 
Taylor expansion 

3 2 r 2 

7(r) ^ 1+ 2^' (6) 

This approximation is valid for not too large values of 3. In this case, Eq. (§) becomes an 
equation for particles moving in a spherical harmonic well with a rigid boundary. 
The stationary states are 

Xe -«2 ; rV2 r * iF A* + 3 / 2 ) ■-„/ , 0/0. ,.2. .2 



1 + 3/2;^], (7) 
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where 1F1 is a confluent hypergeometric function and Yi m / is a spherical harmonic. The 
quantities A n i are the normalization constants and are found from the condition 

<->■ 

J d 3 rr nM (r,t)^mm'(r,t) = l. (8) 
The parameters a n i and e n \ are related to the energy eigenvalues E n i by 
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a A nl = m(E n i - m)(3 2 /R 2 , 

e 2 nl = E nl (E nl - 2m) . (9) 



E ni are given as the solutions to 

i*i + 3 / 2 ' ^ R2 ) = • ( 10 ) 

The normalized contribution to the momentum distribution from the energy state with 
quantum numbers n, I, m' is given in terms of the absolute value squared of the Fourier 
transform of Eq. (|7|) , namely, 

Vw(p) = Jd 3 re-^ nlm/ (r). (11) 

Since the problem has an azimuthal symmetry, the wave function in momentum space does 
not depend on the quantum number m' and is a function only of the momentum magnitude. 

Vw(p) = Vw(p)<Wq • (12) 

Consequently, the momentum distribution is obtained by weighing the contribution from 
each state with the statistical Bose-Einstein factor and adding up the contribution from all 
of the states 

d 3 N _ x , <p nl (p) 

2^ „fK„,- u WT 1 ^ L6 > 



where <fi n i(.p) is defined by 



2E n i 

<i>nl (p) = T^Tg | Sm'Olpnlm' (p)\ 2 (14) 

and the chemical potential /i is computed from 

n, I 

for a given number of particles N. Eq. ( p!5|) follows from Eq. ( pT3| ) after integration over c/ 3 p. 

Fig. 1 shows the systematics obtained by varying the parameters involved. The curves 
in Figs, la, b and c are computed for mid rapidity pions, y = 0, for which the assumption 
of spherical expansion should not be important since these do not experience the effect of 
longitudinal flow. Fig. la shows the behavior of the distribution for a fixed temperature 
T = 120 MeV, a fixed value of the surface expansion velocity f3 = 0.5 and a total number of 
particles iV = 150 for various values of the fireball's radius R. Notice the convex shape of the 
distribution at low m t — m for large values of R and the transition to a concave shape with 
decreasing values of R. This is a density effect since for large R and a fixed total number 
of particles the density is lower than for smaller values of R. In the former case, the value 
taken by the chemical potential is far from the first energy state, whereas in the latter, the 
chemical potential is close to the energy of this state and thus the lowest lying energy states 
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contribute with a more significant statistical weight. The same effect can be obtained by 
keeping a fixed radius and varying the number of particles. 

Fig. lb shows the behavior of the distribution when varying the temperature T main- 
taining fixed values of the radius R = 8 fm, the surface expansion velocity (3 = 0.5 and the 
total number of particles N = 150. As could be expected, the main effect goes into the 
effective slopes describing the distribution's overall fall off. The distribution for T = 100 
MeV rises more steeply at low values of mt — m, due to the proximity of the system to the 
critical temperature for Bose-Einstein condensation, as discussed in Ref. ||. Fig. lc shows 
the effect of varying the surface expansion velocity comparing the cases with (3 = and 
(3 = 0.5 keeping fixed the values of the radius R = 8 fm, the temperature T = 120 MeV 
and the total number of particles N = 150. Notice that the effect corresponds also to as- 
signing different overall effective inverse slopes to each curve, the largest one corresponding 
to the curve with non-vanishing expansion velocity. Finally, Fig. Id shows the shape of the 
distribution for a rapidity away from the central region, in this case with the parameters 
T = 120 MeV, (3 = 0.5, R = 8 fm and iV = 150 for y lab = 3.0. Notice the pronounced 
bending upwards of the distribution at large values of m t — m. This effect can be attributed 
to a larger density of states at high energy eigenvalues as compared to a calculation without 
boundary. 

With these systematics at hand, we proceed to describe the mid rapidity pion data on 
central Au+Au reactions at 11.6 A GeV/c. 0. Data correspond to pions within a rapidity 
interval \Ay\ < 1 around central rapidity. Given our assumption of a spherically symmetric 
fireball, our calculation will compare best to this part of the spectrum since these are the 
pions that do not experience the effects of a longitudinal flow that might be different from 
the flow in the transverse direction. 

Rather than performing an exhaustive search in the whole parameter space and in order 
to test the plausibility of this type of description, here we fix the value of the parameters 
involved to reasonable and more or less accepted values. It is clear that a more complete 
analysis requires the proper treatment of Coulomb effects, this will be the subject of an 
upcoming work. We take T = 120 MeV, f3 = 0.5 (corresponding to an average collective 
expansion velocity < v >~ 0.4c). For the mean negative pion multiplicity we take N n - = 
160 We consider a fireball radius R = 8 fm. Fig. 2a shows the theoretical distribution 
compared to data. In order to compare with the invariant differential cross section reported 
in Ref. which is normalized to a subset of (central rapidity) 116 negative pions, the curve 
has been multiplied by a constant Af = 0.56 that minimizes the x 2 when we compare to 
data above m t — m = 0.4 GeV, the region where Coulomb effects should start becoming 
less significant. Above m t — m = 0.4 GeV the agreement between data and theory is very 
good. Below m t — m = 0.4 GeV the curve follows the shape of the data points but these 
last are still above the calculation. This could be a good feature since one knows that the 
long-range Coulomb effects should push the distribution for low momentum negative pions 
upwards, given that their Coulomb interaction with the overall positive charge is attractive. 

We now use these parameters to describe data on positive pions. Fig. 2b shows the 
theoretical distribution calculated for T = 120 MeV, (3 = 0.5, R = 8 fm but a total positive 



pion multiplicity N n + = 115 [[TJ[], compared to data. In order to compare with the invariant 
differential cross section reported in Ref. which is normalized to a subset of (central 
rapidity) 94 positive pions, the curve has been multiplied by the constant Af = 0.59 that 
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minimizes the x 2 when we compare to data above m t — m = 0.4 GeV. The agreement 
between data and theory is also very good for the region above m t — m = 0.4 GeV. However, 
the raise of the curve below m t — m = 0.4 GeV is less steep than for the negative pion 
case. This is easy to understand since for positive pions the density is lower than for the 
negative ones. Also in the same region, the theoretical curve is marginally below data and 
the Coulomb distortion will push it even more below. We speculate that this signals that 
data prefer a slightly lower value of the radius but again, this can only be confirmed after 
the proper inclusion of Coulomb effects. 

In conclusion, we have shown that a very good description of mid rapidity, charged pion 
spectra can be achieved by a phenomenological calculation whose key ingredient is the proper 
treatment of the large pion density produced in central, relativistic Au+Au reactions. This 
large density leads to consider the pion system as confined during the early stages of the 
collision, before decoupling. Such scheme has consequences on both ends of the spectra. At 
low transverse mass, the convex shape of the distribution is due to the large value of the 
chemical potential associated with the mean pion multiplicity per event. At high transverse 
momentum, the convex shape of the distribution is due to the higher density of states as 
compared to a calculation without boundary. Another important element is the inclusion 
of collective flow. We find that an average flow velocity of < v >~ 0.4c together with 
a temperature T = 120 MeV does a very good job describing data above m t — m = 0.4 
GeV when the fireball's radius is about R = 8 fm. A more conclusive statement can be 
made only after we include a proper treatment of the Coulomb corrections Perhaps 



more importantly is the fact that the different rates at which the positive and negative pion 
distributions rise at low values of m t — m can be understood in part as an effect related to 
the correspondingly different measured yields, making the negative pion subsystem denser 
that the positive pion one. 

Support for this work has been received in part by CONACyT Mexico under grant No. 
I27212-E. 
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FIGURE CAPTIONS 

Fig. 1: Systematics obtained by varying (a) the fireball's radius R, (b) the temperature T 
and (c) the surface expansion velocity (5 for rapidity y cm = 0. (d) Distribution for a rapidity 
Vlab = 3.0. 

Fig. 2: (a) Theoretical distribution {2nm t )~ 1 d 2 N / dm t dy computed with the parameters 
T = 120 MeV, (3 = 0.5, N v - = 160 and y cm = compared to data from the E-802/866 on 
mid rapidity negative pions from central Au+Au reactions at 11.6^4 GeV/c. The theoretical 
curve has been multiplied by the constant Af = 0.56 that minimizes the y 2 when compared 
to data above m t — m = 0.4 GeV. (b) Distribution computed with the same parameters but 
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with N n + = 115 compared to data on mid rapidity positive pions from the same reaction. 
The theoretical curve has been multiplied by the constant M = 0.59 that minimizes the 
X 2 when compared to data above m t — m = 0.4 GeV. Data are {2nm t <Jtri g )~ 1 d 2 a / dm t dy in 
the rapidity interval < Ay < 0.2. The total measured yield spans the rapidity interval 
| Ay | < 1 around central rapidity 0. 
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